Parameterization of the Woods-Saxon Potential for Shell-Model Calculations 
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The drastically expanded use of the Woods-Saxon potential in modern day nuclear physics and the 
availability of new nuclear data motivated us to review and optimize the parameters of this potential 
to the experimental single-nucleon spectra around the doubly-magic nuclei between 16 and 208 Pb. 
We obtain a parameterization which is applicable over the whole nuclear chart for nuclides between 
16 O and the heaviest elements. Apart from Coulomb components the obtained parameter set is 
isospin symmetric. We demonstrate that the potential provides a good description of the nuclear 
mean field leading to quality single-particle spectra, nuclear radii, prediction of drip-lines, shell 
closures and other properties. This presented Woods-Saxon fit provides an adequate single-particle 
basis for shell model calculations bridging over into the continuum. 

PACS numbers: 21.19.-k,21.60.-n 
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I. INTRODUCTION 

The description of the nuclear many-body system by 
an effective mean field is a doorway to the understanding 
of atomic nuclei. The averaged single nucleon dynamics 
in the field of all other nucleons is a starting point in prac- 
tically all many-body methods. In this context, a good 
choice of single-nucleon basis states is the key for suc- 
cess in any quantum many-body approach. A plethora of 
experimental observations such as magic numbers, shell 
gaps, binding energies, nuclear radii, abundances of nu- 
clei in nature and reaction properties all confirm the re- 
markable success of the very simple, pure mean field pic- 
ture. 

Generally, the motion of non-interacting particles in 
the mean field is not an exact solution to the many-body 
problem. Residual interactions or collective dynamics of 
the mean field itself are present. The Hartree-Fock or 
Hartree-Fock-Bogoliubov approach allows to variation- 
ally find the best possible mean field thus minimizing the 
residual interactions. The technique has been demon- 
strated to be very successful. Nevertheless, at present 
it is still a challenge to find the best mean field with 
properly preserved symmetries, not to mention the inter- 
actions and related physics that go beyond the Hartree- 
Fock approach. 

In contrast, in a very simple approach the mean field 
can be taken in the form of a three dimensional har- 
monic oscillator, which provides an analytical set of ba- 
sis states. The possibility of an exact translationally in- 
variant treatment with full center of mass extraction is 
particularly appealing. Historically, the harmonic oscil- 
lator mean field with added spin orbit term was the first 
successful mean field treatment, by which the correct se- 
quence of orbitals and the magic numbers was predicted 
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In order to achieve a quantitative description of nu- 
clei, oscillator parameters have to be selected in depen- 
dence of the nuclear mass. Typically, for spherical nuclei 
of mass A, the oscillator parameter hui = Al/A 1 ^ MeV 
provides a good description, which was found by compar- 
ing the oscillator rms mass radius with experiment. The 
analytical form of the single particle basis comes at the 
expense of relatively strong residual particle-particle in- 
teractions. Still, very successful shell-model descriptions 
have been developed on the basis of the harmonic oscilla- 
tor [2|], where the many-body Hamiltonian that includes 
residual interactions is diagonalized exactly. 

Present day nuclear physics has expanded its reach far 
beyond the valley of stability into regions of the nuclear 
chart where the continuum spectrum of the mean field 
potential becomes important. To address these situa- 
tions or for the purposes of nuclear reaction physics, the 
choice of a Woods-Saxon Q potential as a model for the 
mean field has become a common approach. The ana- 
lyticity of solution in the harmonic oscillator is traded 
for the presence of a continuum spectrum and a poten- 
tial resembling the geometric distribution of nuclear den- 
sity. Although the single particle states must be found 
numerically, such computations are trivial with today's 
computers, making the numerical approach preferential 
over analytically solvable models such as Ginocchio po- 
tential [4( or the square well. Other forms of potentials 
have been considered in the past [f|. With little excep- 
tion, all modern theoretical techniques dealing with the 
physics on the interface of structure and reactions have 
their roots in the Woods-Saxon potential @, 0] • 

A number of parameterizations of the Woods-Saxon 
potential have been published, created with different ob- 
jectives and relevant to different nuclear mass regions. 
Most commonly used is the so-called "Universal" pa- 
rameterization |8j , which was adjusted to reproduce the 
single-particle binding energies of proton and neutron or- 
bitals around the doubly-magic nucleus 208 Pb and correct 
ground state spins for nuclei of masses around A=180, 
but claimed to be applicable to lighter mass regions as 
well. Characteristic for the "Universal" parameterization 
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is the choice of different radii for the proton- and neutron 
potentials. It has been pointed out, that this parameter- 
ization has shortcomings with respect to lighter nuclei 
and that it predicts charge radii inconsistent with exper- 
iment. 0. 

The importance of a good starting parameterization of 
the mean field potential, its expanded use in the context 
of the continuum shell model and the availability of new 
experimental data motivated us to revisit the question of 
parameters in the Woods-Saxon potential. In this work 
we find a set of parameters to be considered as a global 
Woods-Saxon parameterization, obeying the symmetry 
principles of isospin conservation in the nuclear interac- 
tion as well as correct two-body kinematics. We use a 
set of single particle and single hole states in the vicin- 
ity of doubly-magic nuclei 16 0, 40 Ca, 48 Ca, 56 Ni, 100 Sn, 
132 Sn and 208 Pb to provide experimental input on the 
single particle data. We adjust the potential parameters 
to reproduce the experimental data by applying a least- 
squares fit algorithm. 

This work is a continuation and extension of a previous 
effort to determine a globally applicable set of Woods- 
Saxon potential parameters. Among the similarities be- 
tween our work and the previous parameter set is the 
isospin-symmetry of the nuclear potential. While the 
previously obtained parameter set was not published, a 
total-Routhian-surface calculation based on it was dis- 
cussed in jioj . 

The paper is organized as follows. In section[TlJ we de- 
scribe the experimental data that is used in optimizing 
the parameters, in section IIIII we address the properties 
of the Woods-Saxon Hamiltonian, its symmetries and the 
parameterization. In section [TVl we describe the calcula- 
tion of the orbitals of the Woods-Saxon potential and the 
least squares fit algorithm employed to adjust the param- 
eters to fit the experimental data. The predictions of our 
Woods-Saxon potential for rms charge radii and neutron 
radii is discussed in section [Vl as well as its predictions 
for certain neutron-rich nuclides. 



II. SINGLE PARTICLE LEVEL DATA 

Basis for this investigation are the experimental sin- 
gle particle spectra around the doubly magic nuclei 16 O, 
40 Ca, 48 Ca, 56 Ni, 132 Sn and 208 Pb. These systems were 
chosen, because they are expected to show the most 
pure experimental manifestation of single-particle excita- 
tions. The binding energies of the single-nucleon orbitals 
around a magic nucleus of mass A were extracted from 
the excited states of the neighboring nuclei by 



e(A,r) = AM{A)-(AM(A-1)+AM(1))-E(F) (1) 
for hole-states, 

e(A + l,r) = AM(A+1)-(AM(A) + AM(1))+E(r) 

(2) 
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5.288 


4.688 




1.000 



Table I: (£3/2 resonances in O and F, energies are in MeV. 
The bottom line corresponds to the adopted weighted single 
particle energies. 

a The spectroscopic factors are from R-matrix analysis in Ref. 

M 

b Our parametrization with adjusted potential depth was used 
to compute the single particle decay width at the experimen- 
tal energies, the ratio of experimental width to single particle 
determined the spectroscopic factor S. The total strength 
was normalized to 1. 

c The proton channel is only 20% of the total width. 



for particle-states. Here AM(A) denotes the experimen- 
tal mass defects of a nucleus with mass number A and 
E(I V ) the excitation energy of a given state. Both proton 
and neutron orbitals are considered. 

The experimental data were generally extracted from 
the NNDC data base [ll| - [12] , selecting data sets repre- 
senting single-proton or neutron transfer or pickup reac- 
tions, where available. For the 40 Ca, 48 Ca and 208 Pb nu- 
clides, some of the single-particle and single-hole strength 
was observed to be fragmented over a number of individ- 
ual states. In these cases, we calculate the orbital energy 
from an average excitation energy of the states, weighted 
with the single-particle spectroscopic factors. The cases 
for which the data base showed contradictory experimen- 
tal results are discussed separately in the following para- 
graphs. The data set of orbital energies used in this study 
and the relevant experimental information is summarized 
in the Table HO 

In the case of 17 F and 17 0, we used d 5 / 2 and Si/ 2 states 
established by direct transfer-reactions. However, the 
corresponding d 3 / 2 orbits are unbound and their strength 
is fragmented over a number of resonances, both in 17 F 
and 17 O. The resonance data for neutron scattering on 
16 O was taken from Ref. [23j . including the analysis of 
spectroscopic factors. The analogous factors for protons 
were extracted from the ratio of the observed proton- 
resonance widths [IH to a single-particle decay width 
calculated with our Woods-Saxon potential for the ex- 
perimental kinematics. The results are summarized in 
Tab. D 

For the case of proton orbitals in 40 Ca, the Nuclear 
Data Sheets entry for 39 K [l]| shows a state at 8.43 MeV 
excitation energy, which is assigned d 3 / 2 character with 
a spectroscopic factor of C 2 S'=0.24, based on data pub- 
lished in Ref. 24[ . We excluded this state from our anal- 
for two reasons; First it would represent a rather 



ysis, 

unusual fragmentation of single-particle strength, placing 
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a fragment of the ground state at over 8 MeV excitation 
energy. Second, including it would lead to inconsistent 
Coulomb shifts between the proton and neutron orbitals. 
The assignment of d 3 / 2 in Ref. [13] was based on an angu- 
lar distribution measurement, but the difference between 
a assignment and the more likely d 5 / 2 assignment 
seems inconclusive from the data presented in the publi- 
cation. In view of this assessment, we based our analysis 
of the 39 K data on the older measurement of [25] ]. Note 
that, after excluding the above mentioned days assign- 
ment, the data from both references (25| and [24[ leads 
to almost identical energies of s.p. orbitals. 

For 56 Ni, only one experiment for a direct neutron- 
transfer is published, based on the (d,p) reaction [27l ] . 
This experiment revealed spectroscopic factors of single- 
particle character for the lowest 3/2~,5/2~ and 1/2 - 
states. For the proton orbitals, no direct spectroscopic 
factor measurements are available to date, however, 
the lowest 3/2~,5/2~ and l/2~ states show consistent 
Coulomb-shifts, and therefore were assumed to be the 
analog proton excitations. 

100 Sn is still beyond the reach of today's radioactive 
beam facilities, but systematic studies of the shell struc- 
ture of nuclides approaching it have produced a set of 
effective single particle orbitals [30jj . which we include in 
our investigation. 

For the doubly-magic nucleus 132 Sn no single-particle 
reactions were published to date. However, the level 
schemes of nuclei around 132 Sn are well established from 
/3-decay measurements. Using data from the NNDC data 
base [HI [2(3] , we identified the lowest states of given spin 
and parity with the single-particle or single hole states of 
the orbitals in question, and experimental binding ener- 
gies for the ground states. 

208 Pb has been studied by a multitude of experi- 
ments. The spectrum of single-particle and hole excita- 
tions seems well established by the measured large spec- 
troscopic factors [2l|,|22j. Nevertheless, even for this text- 
book doubly magic nucleus, some fragmentation of the 
orbital strengths can be expected for high excitation en- 
ergies. Data on the fragmentation of proton and neutron 
hole states around 2()8 Pb was established in Refs. 28] 
and [2i| . Although these experiments show the main or- 
bital strength to be concentrated in one state for each 
orbital, the smaller fragments lead to significant centroid 
energy shifts, approximately 0.5 MeV for the deepest hole 
states. In contrast to this situation, no data on the frag- 
mentation for the neutron and proton particle states was 
found in the literature, in spite of the large number of 
published experimental studies. It seems worth while to 
revisit experiments on the fragmentation of single par- 
ticle structure around 208 Pb with modern experimental 
techniques in order to address this apparent deficiency of 
the experimental data. 



Table II: Systematics of nuclear single particle energies. In 
the first three columns, we list the orbital type, extracted 
binding energy, and excitation energy of the single particle 
level. The binding and excitation energies are weighted with 
spectroscopic factors (see text), the value of the total spec- 
troscopic strength S = ^ S = C 2 S/(2j + 1) appears in the 
fourth column for the cases where relevant experimental in- 
formation is available. The fifth column identifies individual 
excited states used in the analysis. 



Orbital 


Energy 


Ex 




individual states 




[MeV] 


[MeV] 




[MeV] 


ib 


neutron hole 




ia O[ll|] 




-15.66 










neutron particle 




u 23] 


ld 5 /2 


-4.14 





1.00 





2S1/2 


-3.27 


0.87 


1.00 


0.87 


1^3/2 


1.710 


5.836 


1.00 


see Tab. I 


proton hole 




ib N[n[ 


lPl/2 


-12.13 





1.13 





proton particle 




U F 12] 


ld 5 /2 


-0.60 





0.94 





2Sl/ 2 


-0.11 


0.495 


0.83 


0.495 


W3/2 


4.688 


5.288 


1.00 


see Tab. I 


4U Ca 


neutron hole 




;jy Ca 25] 


W5/2 


-22.39 


6.744 


0.90 


21 states between 4.320 








and 9.500 MeV 


2S1/2 


-18.19 


2.533 


0.82 


2.463, 4.017 


W3/2 


-15.64 





0.94 





neutron particle 




41 Ca 26] 


I/7/2 


-8.36 





0.77 





2P3/2 


-5.84 


2.519 


0.91 


1.940, 2.461, 3.731, 4.602 


2fl/2 


-4.20 


4.157 


0.70 


3.613, 3.943, 4.109, 4.754 


1/5/2 


-1.56 


6.801 


0.95 


24 states between 4.878 








and 9.084 


proton hole 




39 K H| 


ld 5 /2 


-15.07 


6.738 


0.83 


16 states between 5.262 








and 9.750 


2Sl/ 2 


-10.92 


2.593 


0.87 


2.520, 4.095 


W3/2 


-8.33 





0.93 





proton particle 




4i Sc [14] 


I/7/2 


-1.09 





1.12 





2^3/2 


0.69 


1.784 


0.94 


1.716, 2.419 


2P1/2 


2.38 


3.471 


0.75 


3.471 


1/5/2 


4.96 


5.685 


0.33 


5.709, 3.192, 5.862, 6.470 


48 Ca 


neutron hole 




4 'Ca [15J 


ld 5 /2 


-15.61 


5.669 


0.15 


4.980, 5.300, 5.450, 6.250, 








6.870 


2S1/2 


-12.55 


2.600 


0.90 


2.600 


W3/2 


-12.53 


2.580 


0.90 


2.580 


1/7/2 


-10.00 


0.050 


0.85 


0, 3.300, 3.430 


neutron particle 




4H Ca 16] 


2^3/2 


-4.60 


0.545 


0.97 


0, 4.069 


2pi/ 2 


-2.86 


2.282 


1.03 


2.021, 4.261 


1/5/2 


-1.20 


3.946 


0.95 


3.586, 3.993 


lff9/2 


0.130 


5.276 


0.47 


4.018, 4.887, 5.378, 6.529, 








6.753 


proton hole 




47 K [JJ] 
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Table II - Continued 



Orbital 


Energy 


Ex 


S 


individual states 




[MeV] 


[MeV] 




[MeV] 




-21.47 


5.664 


0.62 


3.432, 5.220, 5.465, 6.462, 








7.740,8.020,8.530 


W3/2 


-16.18 


0.377 


1.22 


0.359, 3.930 


2Si/ 2 


-16.10 


0.295 


0.92 


0, 3.850 


proton particle 




4y Sc [a 


l/r/2 


-9.35 


0.278 


0.91 


0, 3.809 


2p 3 /2 


-6.44 


3.187 


0.54 


3.085, 6.717 


2pi/ 2 


-4.64 


4.984 


0.88 


4.495, 5.015, 5.663, 6.816 


bb N j 


neutron hole 




bb Ni [17j 


I/7/2 


-16.64 










neutron particle 




b/ Ni 27] 


2P3/2 


-10.25 





0.91 





1/5/2 


-9.48 


0.768 


0.91 


0.768 


2pi/ 2 


-9.13 


1.113 


0.90 


1.113 


proton hole 




bb Co [nl 


I/7/2 


-7.17 










proton particle 




b 'Cu 18] 


2P3/2 


-0.69 










1/5/2 


0.33 


1.028 




1.028 


2P1/2 


0.41 


1.106 




1.106 


iuu Sn [30], see note" 


neutrons 


2fl/2 


-18.38(20) 








lff9/2 


-17.93(30) 








2t?5/2 


-11.13(20) 








157/2 


-10.93(20) 








3S1/2 


-9.3(5) 








lftll/2 


-8.6(5) 








2d 3 /2 


-9.2(5) 








protons 


I/5/2 


-8.71 








2P3/2 


-6.38 








2pi/ 2 


-3.53(20) 








lff9/2 


-2.92(20) 








2d 5/2 


3.00(80) 








1^7/2 


3.90(15) 









Table II - Continued 



Orbital 


Energy 


Ex 


S 


individual states 




[MeV] 


[MeV] 




[MeV] 




rM S-a, see note " 


neutron hole 


on 


1Q1 


I97/2 


-9.75 


2.434 




9/1^/1 




2d 5/2 


-8.97 


1.655 




1.655 




3Sl/2 


-7.64 


0.331 




0.331 




lftll/2 


-7.54 


0.224 




0.224 




2^3/2 


-7.31 












neutron particle 


133 Sn 




2/7/2 


-2.47 












3P3/2 


-1.57 


0.854 




0.854 




1^9/2 


-0.86 


1.561 




1.561 




2/5/2 


-0.42 


2.004 




2.004 




proton hole 


m In 


19] 


2pi/2 


-16.01 


0.302 




0.302 




159/2 


-15.71 












proton particle 


1M Sb 20] 


I37/2 


-9.68 












2t?5/2 


-8.72 


0.962 




0.962 




2t?3/2 


-6.97 


2.708 




2.708 




1^11/2 


-6.89 


2.793 




2.793 




2U8p b 


neutron hole 




2UVp b 


l/l9/ 2 


-11.40 


4.036 


0.98 


3.400,5.410,5.620 28] 


2/7/2 


-9.81 


2.439 


0.95 


2,340,4.570 |28j 


l«13/2 


-9.24 


1.870 


0.91 


1.630,5.990 28 


3P3/2 


-8.26 


0.89 


0.88 


0.890 


21 


2/5/2 


-7.94 


0.57 


0.60 


0.570 


2M 


3fl/2 


-7.37 





0.90 


21] 




neutron particle 




2U9p b 


M 


2ff9/2 


-3.94 





0.83 







1*11/2 


-3.16 


0.779 


0.86 


0.779 




ljl5/2 


-2.51 


1.424 


0.58 


1.424 




3^5/2 


-2.37 


1.565 


0.98 


1.565 




4si/ 2 


-1.90 


2.033 


0.98 


2.033 




2ff7/2 


-1.44 


2.492 


1.05 


2.492 




3d 3/2 


-1.40 


2.537 


1.09 


2.537 




proton hole 




2U'/rpj 




lff7/2 


-12.00 


3.991 


0.52 


3.469,3.995,4.888 


2d 5/2 


-9.82 


1.806 


0.65 


1.683,4.696 


lftll/2 


-9.36 


1.348 


0.88 


1.348 




2^3/2 


-8.36 


0.351 


0.90 


0.351 




3S1/2 


-8.01 





0.85 







proton particle 




™ M Bi [22| 


I/19/2 


-3.80 





1.00 







2/7/2 


-2.90 


0.897 


1.38 


0.897 




1*13/2 


-2.10 


1.697 


0.93 


1.612,2.601 


2/5/2 


-0.97 


2.824 


0.87 


2.824 




3P3/2 


-0.68 


3.116 


0.98 


3.116 




3fl/2 


-0.16 


3.637 


0.54 


3.637 





a The data was extracted by systematic shell model studies 
of nuclei close to 100 Sn 

b The data is based on the excitation energy of the lowest 
state of given spin and parity, no information on spectroscopic 
factors is available. 
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III. WOODS-SAXON CALCULATIONS 
A. Woods-Saxon Hamiltonian 

In this section, we develop the terms in the Woods- 
Saxon Hamiltonian from very general assumptions about 
the character of the nuclear mean field. For the most 
part, we arrive at the conventional Woods-Saxon Hamil- 
tonian, but the purpose here is to emphasize the aspects 
of the mean field description and to emphasize the the- 
oretical foundations behind the construction. We high- 
light the kinematic aspects of the problem and symmetry 
considerations. The concept of reduced mass and isospin- 
symmetry are of particular importance in the parameter- 
ization of the Woods-Saxon Hamiltonian introduced in 
this work. 

The general strategy behind the construction of a cen- 
ter of mass Hamiltonian for the nuclear mean field starts 
with the assumption of a scalar interaction potential be- 
tween nucleon and core as a sum of the nuclear and 
Coulomb parts. We assume the nucleon and a core form- 
ing a nucleus with mass number A = N + Z containing 
N neutrons and Z protons. Thus the core has A' = A— 1 
nucleons. In the following, we use the prime to denote 
quantum numbers and parameters of the core. 

Woods and Saxon Q suggested to model the nuclear 
mean field i.e. the nucleon-core interaction with a spheri- 
cally symmetric potential that has a Fermi-function form 



f(r,R,a) = 



exp 



r-R 



~] -l 



(3) 



where the size R and diffuseness of the surface a are fixed 
parameters of the same units of length as r. 
The total nuclear potential is defined as 



V(r) = -Vf(r,R,a), 



(4) 



where V represents total strength and the minus sign 
is introduced to represent the attractive nature of the 
interaction. 

The electromagnetic force is a second part contributing 
to the proton-core interaction. This repulsive potential is 
fully determined with the assumption of a given nuclear 
charge distribution p(r). The solution of the correspond- 
ing electrostatics problem gives 

V c (r) = 47re (J- r' 2 p{r')dr' + j r'p{r')dr'^j . (5) 

In the spirit of the Woods-Saxon parameterization it is 
often assumed that the nuclear charge distribution is 
proportional in shape to the same function ([3]) p(r) ~ 
/(r, R c , a c ), where the coefficient of proportionality must 
be determined from the normalization of density to the 
total nuclear charge. The integration in Eq. §5§ along 
with a normalization of density must be done numeri- 
cally, which is often too time consuming. The influence 



of surface terms on the strength of the Coulomb interac- 
tion is, however, weak. We have numerically tested, that 
the diffuseness of the charge distribution can be set to 
zero within the precision of the fit discussed below. Fur- 
thermore, for the same reason we have assumed R c = R 
which removes an extra unnecessary parameter that has 
little influence on the outcome. Except for special cases 
[7| these assumptions are typical in other Woods-Saxon 
parameterizations. Through this paper we adopt the fol- 
lowing form of the Coulomb potential 



V c (r) = Z'e 2 



(3R 2 
1/r, 



! )/(2i? 3 



r <R, 
r>R, 



(6) 



which as a result of the above assumptions corresponds 
to a uniformly charged sphere of radius R, which can be 
treated analytically. 

The understanding of the Woods-Saxon potential as 
a two-body problem naturally leads to the introduction 
of a reduced mass. However, in most previous treat- 
ments, notably [Sj, the bare nucleon masses were used 
in the kinetic energy term of the Hamiltonian. While 
the difference between a reduced or bare nucleon mass 
is not important for the bound state spectrum of heavy 
nuclides, its use would create serious problems for the 
description of light nuclides in our calculations, notably 
the nucleon orbitals around 16 0. Furthermore, we aim 
to apply this potential as a basis for continuum states, 
where the kinematic aspect becomes central. Thus, in 
the kinetic energy operator we use the following reduced 
mass 



A* = 



1 

nil- 



1 

M 7 



(T) 



where is a neutron/proton mass and M' is the mass 
of an A — 1 core, taken to be (A — l)amu 

In addition to the central and Coulomb potentials iden- 
tified above, the relativistic corrections for a nucleon with 
a Fermi momentum of typically 200 MeV are impor- 
tant. In order to highlight the nature of these correc- 
tions, we give a non-relativistic reduction of the two-body 
Coulomb problem to the order of (v/c) 2 . The correction 
to the nuclear potential is discussed below. The Dirac 
Hamiltonian reduces to 

ff = |!-|l +W + 5 L..[V^(r)xp] + i Lv^(r). 

(8) 

Throughout this section, we use natural (Planck) units 
where h = c = 1. In nuclei the contribution from the 
the second (kinetic energy) p term is small and there- 
fore ignored. It turns out that, in comparison to nuclear 
forces, corrections due to Coulomb given by the last two 
terms are also small. For example the last so-called Dar- 
winian term is constant, and only non-zero inside the 
uniformly charged sphere. For a typical nucleus its value 
-3Z'e 2 /(8R 3 p 2 ) is only of the order of 20 - 30 keV which 
should be compared with ~ 50 MeV of a typical total 
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depth V in ([4]). It should also be stressed that terms 
of this nature can be subsumed in the phenomenologi- 
cally determined parameters of the central potential. The 
spin-orbit contributions due to the Coulomb potential in 
eq. [8] are an order of magnitude smaller than the nuclear 
contributions and will be ignored. 

Based on the previous discussion, only the nuclear 
spin-orbit contribution has to be included in descriptions 
of effective nuclear forces in addition to V c . Eq. |(SJ) is not 
valid as a relativistic reduction of the nuclear mean field 
problem due to the complicated structure of the nucleon- 
nucleon force. However, symmetry considerations sug- 
gest the same general proportionality of the spin orbit 
force to be the gradient of the mean field potential. Thus, 
the total effective Hamiltonian becomes 

H=£ + V(r) + V e {r) + ^(jLv{r)) (9) 

where - unlike for the Coulomb field - the potential V{r) 
is not equal to the original potential V(r) and may have 
a different form factor [3l|. Therefore, the form factor of 
V(r) is another assumption that goes into construction 
of the Woods-Saxon Hamiltonian 

V(r) = Vf(r,R so ,a S o)- (10) 

Here, Rso an d aso stand for the radius and the diffuse- 
ness of the spin-orbit term. 

In principle all symmetry preserving forces which in- 
volve the single particle operators p, r, s, t and the core 
spin and isospin operators T' and I' can appear in the 
particle-core parameterizations. These generally small 
terms are extensively discussed and studied within the 
optical model approach to reaction physics [H, [33[ . The 
consideration of an odd-particle mean field that carries 
quantum numbers of an unpaired nucleon is beyond the 
goals of the parameterization discussed here. Therefore, 
we assume that the core generally carries no spin degree 
of freedom I' = 0. The only time such terms can be of rel- 
evance is when single-particle energies are to be extracted 
from hole states Eq. |T]) in this case the particle spin may 
couple to identical single-hole quantum numbers of the 
core. These terms are generally small and have a 1/A 
dependence @, just as in Eq. (TP3"|) . The only "second- 
order" term of concern is a possible isospin-dependence 
of the spin-orbit strength i.e. (T • t)(l • s) This depen- 
dence has been the subject of several studies [HI HU as 
it represents one of the most interesting, fundamental 
and at the same time controversial questions related to 
the nuclear mean field. The example of a non-relativistic 
reduction of the Coulomb potential and the more rigor- 
ous studies using the Relativistic Mean Field approach 
suggest the same sign of isovector term in spin orbit po- 
tential as in the central. However, a substantial amount 
of experimental data [H, [35|, [36| as well as theoretical 
understanding of an issue based on the Walecka model 
suggest the opposite. Namely, for the neutron rich nuclei 
the spin orbit splitting for the neutrons is larger than for 



the protons. Those studies suggest the value kso ~ —0.3 
which can be compared to the typical strength of the iso- 
topic term in the central potential k ~ 0.6 to 0.9; see 
discussion below for details and definitions. The more 
in-depth arguments and understanding of the spin orbit 
splitting comes from studies of shell evolution due to the 
tensor forces [53|- The effects of the p meson exchanges 
in the nucleon-nucleon interaction are particularly im- 
portant for these mean field properties. 

We used our data set and fit procedure to assess the 
variation of spin-orbit strength with isospin. Unfortu- 
nately, as was already pointed out in Ref. [H, [36| 
the data based on single-particle energies is only weakly 
sensitive to this parameter. Within the precision of 
our study we were unable to distinguish any significant 
isospin dependence of the spin orbit term. Thus, with no 
better alternative we define the spin orbit strength in our 
parameterization with no explicit dependence on isospin. 



B. Woods-Saxon parameterization 

The Hamiltonian of the model is defined in Eq. ^j, 
where the central potential is determined in Eq. (TJJ and 
the spin-orbit term is given by Eq. (jTUJ) . Both central 
and spin-orbit parts have a Woods-Saxon form factor ([3]). 
The Coulomb part is given by Eq. ^ . For each individ- 
ual nucleus the Hamiltonian is determined by the fol- 
lowing list of parameters R, a, V, Rso, a so, V and the re- 
duced mass /i; these seven parameters define the Woods- 
Saxon potential. 

The parameters of the potential change as one goes 
over the nuclear chart. The dependence of the above 7 
parameters on the number of protons and neutrons in the 
core-nucleon problem defines the structure of parameter- 
ization. In the following, we will describe our choice of 
structure parameterization, which we call the "Seminole" 
parameterization. 

In the conventional parameterization of the Woods- 
Saxon potential, as well as in our work, the size of nuclear 
potential is calculated as 

R = R C = RoA 1/3 , Rso = Ro^soA 1/3 , (11) 

in terms of the parameters Rq and Rq.so, which are con- 
stant over the nuclear chart. 

The surface diffuseness for both the central and spin- 
orbit potential is assumed to be constant and size inde- 
pendent 

a = a so = const. (12) 

While the isospin dependence of the Coulomb force is 
obvious, the behavior of the effective nuclear potential 
on the isospin of the nucleon t and the core T' has to be 
introduced phenomenologically. We adopted the sugges- 
tion by Lane [37j ] to introduce an isospin dependence to 
the potential by the lowest order isospin invariant term 
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K = V (l-^<t-T'>). (13) 

We chose the "minus" sign so that the parameter k will 
be consistent with conventions. For the ground-state of 
a nucleus, the isospin quantum number is T — \T Z \ — 
\N - Z\/2, which together with the relation t + T' = T 
leads to 

( 3 N = Z 

-4(t ■!") = < ±(N - Z +l) + 2 N>Z , (14) 
[ ± {N-Z-1) +2 N < Z 

where here and below we use upper sign for the proton 
and the lower sign for the neutron. Traditionally, the 
isospin-dependence of the Woods-Saxon potential had 
been parameterized by the expression. 

v = v^ ± .^l) (15 ) 

For heavy nuclei with large neutron excess, the difference 
between the two definitions is small. However, the defini- 
tion of Eq. [14] leads to significantly different predictions 
in lighter nuclides around N=Z. We discuss the validity of 
this assumption for the description of our data in section 

El 

The structure of the "Seminole" parameterization is 
given by Eqs. (TIB")) and (fT4")) which define the dependence 
of the potential depth on Z and iV in an isospin conserv- 
ing way. 

The spin-orbit interaction strength is determined by 

V = AVb, (16) 

using a proportionality constant A which is a constant, 
making V a constant. This aspect of our parameteri- 
zation is different from the traditional approach, where 
the spin orbit potential depth is modified with the same 
isospin dependent factor as the central potential. Finally, 
the reduced mass is given as in Eq. l[7|. where the mass 
of the core is parameterized as (A — l)u. The 6 constants 
Vq, i?0j Ro,SO, a — a 50 5 A, and k and are the actual pa- 
rameters of the global "Seminole" Woods-Saxon fit. 

We have also studied a possible variation of the poten- 
tial radius with isospin. From our investigation of nuclear 
charge radii, which will be described in section TlV C[ we 
concluded that a modification of potential radii with re- 
spect to the standard structure of parameterization in 
Eq. (|lip would lead to wrong predictions of experimen- 
tal charge and neutron radii. 

C. Existing parameterizations 

A number of parameterizations and parameter sets are 
available in the literature. We list the most commonly 
used ones in Table [TTTJ All of these parameterizations use 



Eq. (|15p for the isospin dependence of the nuclear poten- 
tial. In contrast to the "Seminole" parameterization, the 
same isospin dependence enters into spin-orbit term via 

V = XV, (17) 

which should be compared with Eq. (TT()|) . 

The "Rost" parameters [38| were determined from the 
orbital energies of 208 Pb. The "Optimized" [H pa- 
rameter set took the central potential parameters from 
"Rost" , but changed the spin-orbit interaction in order 
to improve predictions for high-spin spectra in the lead 
region. A further refinement of these parameters was 
introduced as the "Universal" parameter set, which im- 
proved the description of high spin states in 146 Gd Q. 

Common to these three parameterizations is the as- 
sumption of different parameters sets for the proton and 
neutron spectra, and specifically a larger neutron than 
proton potential radius. The implementation of these 
parameterizations through the computer code SWBETA 
[45j | uses the the bare nucleon mass in the kinetic energy 
operator and a reduced mass \x — u in the spin-orbit 
coupling potential. 

IV. OPTIMIZATION OF WOODS-SAXON 
POTENTIAL PARAMETERS 

A. Calculation of the Woods-Saxon spectrum and 
adjustment of parameters 

Special attention had to be paid to the way in which 
the experimental particle and hole states are compared to 
the Woods-Saxon calculation. Single neutron states are 
naturally calculated as the potential for (N m + l,Z m ), 
and single proton states as the (N m ,Z m + 1), where 
N m and Z m are the neutron and proton numbers of 
the respective magic core. The hole states around the 
same magic nucleus are nucleons moving in the field of a 
core, which has one nucleon less than the magic core, so 
that the orbitals based on hole states are calculated at 
{Nmi Z, n ). 

As a primary computational tool, we used a Woods- 
Saxon code which is a part of the Continuum Shell Model 
code CoSMo discussed in [42j and that has been success- 
fully used in studies of realistic nuclei [43|. Web imple- 
mentation of this code for bound states along with the 
Seminole parameterization obtained as a result of this 
work can be found in [44j. The spectrum, radii, struc- 
ture of wave- functions and reaction calculations were per- 
formed using this code which makes the current Woods- 
Saxon parameterization readily available for further in- 
corporation into the Continuum Shell Model. 

The independent use of the publicly available program 
SWBETA [45| to calculate the spectrum of the Woods- 
Saxon Hamiltonian provided a crucial verification. This 
program was modified slightly to incorporate the changes 
to the Woods-Saxon parameterization we introduced. 
Note that, while SWBETA calculates the spectrum of a 
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Parameterization 




V [MeV] 


K 


Ro [fm] 


a [fm] 


A 


Roso [fm] 


'Rost'[38j 


n 


49.6 


0.86 


1.347 


0.7 


31.5 


1.28 




P 






1.275 




17.8 


0.932 


'Optimized' [39] 


n 


49.6 


0.86 


1.347 


0.7 


36 


1.30 




P 






1.275 




36 


1.30 


'Universal' [8] 


u 


49.6 


0.86 


1.347 


0.7 


35 


1.31 




P 






1.275 




36 


1.32 


'Chepurnov'[40] 




53.3 


0.63 


1.24 


0.63 


23.8 * 


1.24 


'Wahlborn'[41J 




51 


0.67 


1.27 


0.67 


32 


1.27 



Table III: Commonly used Woods-Saxon parameter sets in the literature. * The "Chepurnov" parameterization introduced an 
additional isospin-dependence for the spin-orbit interaction in deviation from the commonly used proportionality to the central 
potential. 



deformed Woods-Saxon potential, zero deformation was 
used for all calculations presented here. 

We also developed a parameter fit program, which calls 
the Woods-Saxon calculation through a shell command 
with the hypothetical fit parameters and reads back the 
calculated orbital spectrum. The actual parameter opti- 
mization is realized through a Levenberg-Marquardt al- 
gorithm [4t| , which minimizes the x 2 between the exper- 
imental and calculated orbital energies by systematically 
varying the potential parameters. 

While the spectrum of all nuclides listed in section |TT] 
was usually fit with the same set of parameters, the fit 
program also can allow for certain parameters to vary 
for each individual nuclide. This option was used to sys- 
tematically investigate how to improve the dependence of 
certain parameters on the nuclide mass or the neutron- 
proton asymmetry. 

Armed with a fitting procedure we first consider a 
structure of the "Universal" parameterization and fit its 
nine parameters to our data set of orbital energies around 
16 0, 40 Ca, 48 Ca, 56 Ni, 100 Sn, 132 Sn and 208 Pb. The re- 
sulting parameters called "Fit A" are listed in table IIVI 
together with results from the "Universal" parameters. 
We observe that the large difference in central potential 
radius between neutron and proton radii disappears in 
the fit, although these parameters were free to adjust in- 
dependently. The spin-orbit strength parameter A is the 
only parameter in "Fit A" displaying substantially asym- 
metric values for neutrons and protons. Table IIVI also 
lists the RMS energy deviation obtained for the different 
nuclides. The "Fit A" parameters show a more consis- 
tent description of the lighter nuclei than the "Universal" 
parameters and a better description of the nuclides with 
N=Z. The features of "Fit A", namely a symmetric pa- 
rameterization for protons and neutrons, serve as a guide 
to find an improved parameterization and parameter set 
addressed below. "Fit B" , which is included in the table 
will be discussed in the following section. 



B. Changes to the parameterization 

In this paragraph we investigate the validity of 
the isospin-conserving parameterization introduced in 



Eq. HU). In the procedure listed as "Fit B" in table 
H*Vl we introduced the reduced mass of Eq. (fTJ), and the 
constant spin-orbit potential strength of equation (|16p . 
The central potential depth was optimized with a sep- 
arate value for each neutron and proton, particle and 
hole data set in each nuclide, a total of 28 values. The 
other parameters Rg, Ro,so a and aso were free to ad- 
just one common value for all nuclides. The extracted 
V values are displayed in FigO] plotted as a function 
of the isospin multiplier —4 (t ■ T'). The values show 
a linear dependence on the isospin factor, in agreement 
with the parameterization (|14[) . The only significant de- 
viations from the linear dependence are found for both 
the particle (| (t • T'} = 0,V = 49.8 MeV ) and hole 
orbitals (| (t • T'} = 0.1875, V = 57.3 MeV) of Oxygen. 
These data points are an indication that for lighter nuclei 
a reduced potential depth would improve the fit quality. 
However, we did not introduce a mass dependence to the 
potential depth parameter for the sake of simplicity and 
a more stable extrapolation. 



> 




-4/A < t ■ T > 

Figure 1: Central potential strength extracted from "Fit B" 
as a function of the isospin factor - A (t-T'} (see Eq.[l4]|. The 
straight line represents the central potential parameterization 
assumed in this work. 

The question, whether the spin-orbit potential depends 
on isospin was investigated by allowing a modification 
in analogy to the isospin-dependence of the central po- 
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Vb[MeV] 


K 


Ro [fm] 


a [fm] 


A 


Roso[im] 


16 Q 


RMS enerj 
40 Ca 48 Ca 


U deviation [MeV 
56 Ni 100 Sn 132 Sn 


] 

208p b 


Univ. 


u 


49.6 


0.86 


1.347 


0.7 


35 


1.31 


2.65 


1.68 


1.07 


1.18 


1.20 


0.40 


0.37 




P 






1.275 




36 


1.32 


2.31 


2.18 


1.63 


1.40 


2.29 


0.61 


0.43 


Fit A 


n 


52.78 


0.636 


1.282 


0.62 


24.0 


1.18 


1.91 


0.84 


1.06 


0.55 


0.72 


0.53 


0.60 




P 






1.288 




21.2 


1.05 


1.96 


0.69 


0.64 


0.44 


1.31 


0.43 


0.38 


Fit B 




var 


0.000 


1.260 


0.68 


25.5 


1.21 


0.33 


0.52 


1.34 


0.45 


0.74 


0.53 


0.40 






var 












0.32 


0.58 


0.76 


0.39 


1.19 


0.34 


0.32 



Table IV: Parameters and fit quality of the "Universal" parameterization compared to parameter fit "Fit A" adjusting the same 
parameters to the orbital energy data presented in section[TU In "Fit B" parameters labeled "var" were allowed individual values 
for each nuclide for protons, neutrons, for particle and hole state (see text). The individual parameter values are displayed in 
Fig. [T] The values of the other parameters were allowed to adjust one common value during the fit, with the exception of k, 
which was held constant at 0. 



tential. Note that in the previous parameterizations [8(, 
the spin orbit potential strength was assumed to be di- 
rectly proportional to the central potential, thus having 
an explicit isospin-dependence. As was stated earlier the 
isospin-dependence of a spin-orbit term, although an in- 
teresting question, can not be fully addressed with our 
approach. All our attempts to adjust this parameter led 
to a slight deterioration of the fit with most effects due to 
the (T-t)(l-s) term being subsumed in the adjustment of 
the other parameters. Therefore, we defined our param- 
eterization as not containing this term. We believe that 
this feature given its marginal and elusive effect on the 
single particle bound and resonant spectra is not essen- 
tial for our purpose. In fact the isospin dependence in the 
spin orbit is a weak manifestation of the physics that goes 
beyond the mean field and thus it should be treated using 
a more advanced approach such as shell model based on 
the Woods-Saxon basis from our parameterization. Fur- 
thermore, the reduced number of free parameters allows 
us to achieve a more stable parameterization. 



C. Charge Radii 

While the description of the orbital energy spectra is 
an important component of any model of single particle 
motion, the Woods Saxon potential also strives to re- 
produce the geometry of the nucleus and properties that 
are sensitive to the geometry of its wave functions. Es- 
pecially since the radial and potential depth values pa- 
rameters show strong correlations in their influence on 
the particle energy spectrum, it is important to study 
the geometric properties of the wave-functions in an in- 
dependent fashion. 

We took data on experimental rms charge radii from 
the compilation by I. Angeli |47|, for the proton-magic 
nuclides 16 - ls O, 40 Ca - 50 Ca, 108 Sn - 115 Sn and 190 Pb 
- 213 Pb. Using our new parameterization, we calculated 
the RMS radii of the same nuclides by averaging in square 
all proton orbital radii with their proper occupation num- 
bers up to the Fermi level. To the orbital RMS radius, 
the RMS radius of the proton was added in square by 

Rcharge = \l ^orbit ^proton- 



The experimental charge radii, the results obtained 
for our parameters and the results for the "Universal" 
parameters are displayed in Fig. [5J Note that for each 
isotopic chain displayed here, the protons occupy the 
same orbitals, so that the observed variations are due 
to changes in the orbital geometry only. 

We varied the potential radius parameter for values be- 
tween 1.25 and 1.29 fm, while adjusting the other poten- 
tial parameters for each step to reproduce the orbital en- 
ergy data. The value of Rq = 1.26 fm best reproduces the 
charge radii for the 190 Pb- 208 Pb isotopes. The increase 
of experimental charge radii observed beyond 208 Pb has 
been discussed in the context of various models, most 
recently in a systematic discussion of quadrupole cor- 
relation effects in a self-consistent calculation with the 
SLy4 interaction (48j . Our much simpler model can not 
reproduce this increase, either. The Tin and Calcium 
isotopes exhibit slightly larger rms radii than calculated. 
It is likely that this deviation is due to an onset of defor- 
mation in mid-shell, consistent with the calculations of 
0. 

It should be noted, that the variation of rms radii 
within the isotopic chains is determined only in part 
by the global scaling of the potential radius with A 1 / 3 . 
Clearly evident in the data is a systematic dependence 
on the binding energy of the fermi orbital, where proton- 
rich nuclides have relatively larger charge radii than ex- 
pected from the global scaling of the potential radius. 
Our calculations illustrate, that the same orbitals occupy 
a larger fraction of the potential radius as they become 
less bound. 

Direct data on the neutron radii is much harder to 
obtain; However, some information on the relative radii 
of neutrons and protons can be deduced from isovector 
spin-dipole resonances, studied e.g. with the ( 3 He, t) re- 
action at intermediate energies. In Ref. [4!|, the authors 
extracted values for the difference between neutron and 
proton rms radii for even-mass Tin nuclides, which we 
display in Fig. [3J along with the values calculated for 
our Woods-Saxon parameterization and the "Universal" 
parameterization. The data presented in the figure sup- 
ports the assumption of symmetric proton and neutron 
potential radii rather than the relatively larger neutron 
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Figure 2: Experimental rms charge radii divided by the scaling factor A 1//3 as a function of mass number (symbols), for a set 
of Oxygen, Calcium, Tin and Lead isotopes. The data is compared to the results calculated with our parameters (solid line) 
and the "Universal" parameters (dashed line). The experimental uncertainties are typically smaller than 0.005 fm. 



potential of the "Universal" parameterization. 

V. THE NEW WOODS-SAXON POTENTIAL 
PARAMETERIZATION 

A. The Parameters 

We arrive at the structure parameterization given by 
Eqs. (O, (fTTj) - (fl"6|) . The corresponding parameters were 
adjusted to the orbital energy data. Their values are 
displayed in table IVl 



V [MeV] 


K 


Ro [fm] 


a = aso [fm] 


A 


Ro, so [fm] 




52.06 


0.639 


1.260 


0.662 


24.1 


1.16 





Table V: Woods-Saxon potential parameter set obtained in 
this work. 



Summarizing, the important properties of the param- 
eterization and the parameter set are: 

• Reduced mass in the form given by Eq. |[7J) as- 
sures correct two-body kinematics of nucleon and 



core. The average mass per nucleon in a nucleus 
is assumed to be equal to an atomic mass unit 
u = 931.49 MeV/c 2 . The introduction of the re- 
duced mass led to a substantial improvement in the 
quality of description for the Oxygen and Calcium 
spectra. 

• Identical potential parameters for neutrons and 
protons, which in conjunction with Eq. (|13[) assures 
isospin conservation. 

• The spin-orbit term is constant and does not de- 
pend on isospin. 

• There are only six parameters that were fitted. 
These parameters are given in table IVl 

B. Discussion of Fit Results 

The experimental and calculated orbital energy spec- 
tra are displayed in figures l4lfl0l We will first address the 
spectra of the symmetric N=Z nuclides 16 0, 40 Ca, 56 Ni 
and 100 Sn. For these nuclides, the spectrum calculated 
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Figure 3: Values for the difference between the neutron and 
the proton rms radii in Sin isotopes, extracted from exper- 
imental data in spin-dipole resonances [49( ] (symbols). The 
data is compared to the results calculated with our fit pa- 
rameters (solid line) and the "Universal" parameters (dashed 
line) . 



with the "Universal" parameterization shows consistently 
too small neutron and proton energy gaps around the 
magic numbers 8,20,28 and 50. The new parameteriza- 
tion greatly improves the accuracy of the orbital spec- 
trum description. 

As mentioned in section IIIIl the "Universal" param- 
eterization was primarily developed from the spectrum 
of 208 Pb. It is therefore not surprising that it is less 
applicable to symmetric N=Z nuclides, especially since 
it is based on substantially different neutron and pro- 
ton potential parameters. The "Universal" parameters 
fare better with the the neutron-rich nuclides 48 Ca, 132 Sn 
and 208 Pb, which all have similar neutron-excess values 
of (N - Z)/A = 0.167,0.242 and 0.212, respectively. In 
general, the description of these nuclides by the new pa- 
rameters is of similar quality as the "Universal" parame- 
terization, with one exception - although the RMS devi- 
ation of the neutron particle spectrum of 208 Pb with the 
new parameters (0.5 MeV) is similar to the one obtained 
from the "Universal" parameters (0.37 MeV), the new 
parameter spectrum calculates all neutron "particle" or- 
bitals systematically less bound than the experimental 
values. The magnitude of these deviations is not beyond 
the typical discrepancies of individual orbital fits in other 
nuclides, so that we do not think it represents a system- 
atic problem within our parameterization. It should also 
be noted, that we could not find experimental informa- 
tion on the fragmentation of the single-particle strength 
for the nuclide in question, 209 Pb, in the literature. If 
we assume that the orbital fragmentation of high-lying 
states in 209 Pb is similar to that of the neutron hole or- 
bitals in 207 Pb, the systematic "under-binding" of our 
description would vanish. 



VI. PREDICTIONS OF SHELL STRUCTURE IN 
EXOTIC NUCLIDES 

The discussion of all consequences of our new param- 
eterization for the Woods-Saxon potential is clearly be- 
yond the scope of this paper. The reader may further 
explore the features of this parameterization using our 
web-based program [44J. Here, we want to briefly discuss 
the general properties of the shell structure for exotic 
nuclides and the limits of nuclear binding it implies. We 
use the "Seminole" parameterization to calculate shell 
gap energies by filling the orbitals to the Fermi level and 
calculating the ernergy gap to the next unoccupied or- 
bital for all nuclides between Z=8 and Z=150. Table EJ 
lists the nuclides, for which a shell gap of at least 1.5 MeV 
is present for both neutrons and protons and for which 
both protons and neutrons were calculated to be bound. 
Note, that in this table, the gap is calculated from the or- 
bital energy differences within one potential calculation, 
while the data represented in each of the figures HHTUl 
corresponds to two calculations, one for the core and one 
for the odd-nucleon neighbor, which leads to slightly dif- 
ferent shell-gap values. Table |VI] also lists the neutron 
and proton separation energy, calculated as the binding 
energy of the Fermi level. Experimental values are taken 
from the compilation of 50]. 

The data for nuclides between Z=8 and Z=51 is also 
represented as a nuclear chart in Figure [TT1 For nuclides 
with doubly closed sub-shells, the minimum of the re- 
spective proton and neutron shell gaps is represented by 
shades of red. The nucleon numbers of sub-shell clo- 
sures are indicated by thin lines. All other nuclides are 
represented by shades of cyan color, which indicates the 
minimum of separation energy in protons and neutrons. 
We also draw the calculated limits of nuclear binding. 

We find it necessary to stress again that Woods-Saxon 
parameterization is a very crude approximation to the 
nuclear mean field; The true potential can be quite dif- 
ferent (the possibility of deformation is particulary note- 
worthy) and could include more complicated spin and 
tensor structures. Furthermore the parameters of poten- 
tial can change as a function of N and Z in a non- uniform 
manner [5lj |. Although our potential is intended for use 
as a part of a more involved many-body approach it is 
still instructive to explore the shell evolution due to mere 
geometric modification that appear in our parameteriza- 
tion. 



A. N=16 and N=20 Shell Closures 

The modification of shell structure in exotic nuclides is 
currently a topic on which great experimental and theo- 
retical efforts are expended. Otsuka and co-workers have 
attributed a significant renormalization of shell model 
single-particle energies to scalar [52j and tensor [53[ com- 
ponents of the spin-spin proton-neutron residual inter- 
action. These interactions act mainly on the spin-orbit 
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Figure 11: Nuclear chart representing the calculated shell gaps, nucleon separation energies and drip-lines of nuclides between 
nuclides of Z=8 and Z=51. For doubly sub-magic nuclei, which are located at the crossing points of the dotted lines, the 
minimum of the proton and neutron shell gaps is represented by the shade of red, where larger gaps correspond to a deeper 
color. All other bound nuclides are represented by a shade of cyan color, which represents the minimum of the proton and 
neutron separation energies. The calculated proton and neutron drip-lines are also drawn. 
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WS 


exp 


exp 




N 


Z 


Gapjv 


Gapp 


Sjv 


Sp 


Sjv 


Sp 








[MeV] 


[MeV] 


[MeV] 


[MeV] 


[MeV] 


[MeV] 




6 


8 


5.64 


7.01* 


23.80 


5.57 


23.18 


4.63 


16q 


8 


8 


8.97 


8.56 


16.93 


13.24 


15.66 


12.13 


22Q 


14 


8 


2.23 


9.58 


5.58 


24.78 


6.85 


23.26 


24q 


16 


8 


3.51* 


9.68 


3.29 


27.78 


3.6 


26.6** 


28Q 


20 


8 


3.82* 


9.75 


0.23 


32.88 






28 Si 


14 


14 


3.81 


3.71 


16.22 


10.67 


17.18 


11.58 


30 Si 


16 


14 


2.39 


3.95 


10.58 


13.27 


10.61 


13.51 


34 Si 


20 


14 


5.12 


4.24 


7.94 


17.94 


7.53 


18.72 


42 Si 


28 


14 


2.23 


4.44 


3.20 


25.50 


3.20** 


24.56** 


30g 


14 


16 


3.97 


2.46 


18.78 


4.51 


18.98 


4.40 


32g 


16 


16 


2.04 


2.02 


14.23 


8.07 


15.04 


8.86 


40 Ca 


20 


20 


5.95 


5.43 


16.15 


8.66 


15.64 


8.33 


48 Ca 


28 


20 


3.53 


5.95 


8.71 


16.39 


9.94 


15.86 


52 Ca 


32 


20 


1.67 


6.10 


5.02 


19.67 


4.72 


17.80** 


60 Ca 


40 


20 


2.96* 


5.95 


2.53 


24.96 






™Ca 


50 


20 


2.22* 


5.63 


0.25 


30.25 






56 Ni 


28 


28 


4.44 


4.35 


16.01 


6.79 


16.64 


7.17 


68 Ni 


40 


28 


3.29 


4.97 


8.16 


14.86 


7.79 


15.69 


78 Ni 


50 


28 


3.55 


5.13 


5.04 


20.30 


5.62** 




88 Sr 


50 


38 


4.53 


1.57 


10.50 


10.62 


11.12 


10.61 


80 Zr 


40 


40 


2.92 


2.29 


15.49 


3.20 


16.23 


4.44 


90 Zr 


50 


40 


4.66 


2.06 


11.53 


8.19 


11.97 


8.35 


96 Zr 


56 


40 


1.57 


1.96 


6.68 


10.91 


7.86 


11.52 


122 Zr 


82 


40 


3.98* 


1.69 


3.38 


20.60 






100 Sn 


50 


50 


5.15 


5.61* 


16.74 


2.87 


17.65** 


2.80* 


132 Sn 


82 


50 


4.75 


5.23 


7.19 


15.16 


7.31 


15.71 


146 Gd 


82 


64 


5.27 


1.69 


12.04 


4.72 


11.2 


5.38 


182p b 


100 


82 


1.53 


4.43* 


11.04 


2.01 


11.75 


1.31 


208p b 


126 


82 


3.88 


3.64 


7.69 


8.00 


7.37 


8.01 


218 u 


126 


92 


3.82 


1.16 


9.81 


1.78 


8.85 


2.43 


256-fj 


164 


92 


1.37 


1.68 


5.06 


9.64 






276 u 


184 


92 


2.39 


1.81 


2.68 


13.11 






278 114 


164 


114 


1.98 


1.97 


9.28 


2.19 






298 114 


184 


114 


2.84 


1.87 


6.10 


5.41 






342 114 


228 


114 


2.09 


1.70 


3.14 


11.50 







Table VI: Table of nuclides with a Fermi level shell gap calcu- 
lated larger than 1.5 MeV in both protons and neutrons by the 
"Seminole" parameterization. Cases where the upper orbital 
for the shell gap was calculated to be unbound are marked 
with an asterisk. The following columns list the calculated 
neutron and proton separation energies and the experimen- 
tal separation energies from Values marked with two 
asterisks are based on extrapolation, not experiment. More 
calculations may be performed using our program at [44| . 



partners of protons and neutrons. For instance, as the 
proton d 5 / 2 occupation goes from 6 protons at i4Sii6 to 
zero for the g 4 0i6, the effective orbital energy for the 
neutron d 3 / 2 orbital is calculated to increase. This effect 
creates a N=16 shell gap at the neutron drip line and 
also leads to the breaking of the N=20 magic number in 
the "island of inversion" around 32 Mg. 

Our model naturally includes halo phenomena of 
weakly bound states, as observed in the charge radii dis- 
played in Fig. [21 In consequence, the increased radial 
extension leads to modifications in the spin-orbit split- 



ting. The corresponding (sub)-shell gaps for the above 
mentioned example are included in table I VII - the cal- 
culated N=16 shell gap increases from 2.43 MeV in 30 Si 
to 4.51 MeV in 24 0, while the N=20 gap decreases from 
5.09 MeV to 3.82 MeV between 34 Si and 28 0. 

We conclude that some modifications of shell struc- 
ture in exotic nuclides can be explained by a model as 
simple as the presented Woods-Saxon potential, as it in- 
cludes quantum-mechanical dynamics of weakly bound 
particles. In our results however, the size of the N=16 
shell gap does not reach the ~ 6 MeV calculated in . 



B. N=28 and Z=14 Shell Closures 

The question, whether also the N=28 shell closure is 
broken in neutron-rich nuclides is another topic of intense 
investigations and discussion. Here, the mechanics of a 
possible breaking of the N=28 gap would be different 
from the N=20 case discussed above - since the N=28 
closure is based on the orbital and the proton / 5 / 2 
remains empty for the neutron-rich as well as the stable 
N=28 nuclides. 

The Woods-Saxon potential calculation shows a clear 
decrease of the N=28 gap as we compare the values 
for the increasingly exotic 56 Ni, 48 Ca and 42 Si, at 4.44 
MeV, 3.53 MeV and 2.22 MeV, respectively (see Tab ED). 
Based on expectations of a decreasing N=28 shell gap, 
several theoretical investigations [54|- 57] had predicted 
a strongly deformed ground state in f 2 Si 2 s. However, in 
an experiment of one- and two-proton knockout reactions 
with a 44 S beam, it was shown that there is a substantial 
Z=14 shell closure (58l . [59[ and nearly degenerate S1/2 
and g?3/2 orbitals. The shell gap data listed in table IVT1 
and represented in Fig. [11] supports these findings; The 
calculated Z=14 sub-shell gap increases from 3.68 MeV 
in the stable 28 Si to 4.42 MeV in 42 Si. Our calculation 
also shows the near degeneracy of the Si/ 2 and d 3 / 2 or- 
bitals, which here are calculated only 0.63 MeV apart. 
Note that this orbital spacing is also present in the 48 Ca 
data displayed in figure [5] This spacing generates a large 
Z = 14 gap and a very small Z — 16 gap, which is con- 
sistent with a magic character of 42 Si with a stabilized 
N=28 closure and a moderate deformation 44 S, as ob- 
served in [6(J [fil| ■ 



C. Limits of Nuclear Binding 

The limits of neutron binding and location of the neu- 
tron drip-line is a very important property of atomic nu- 
clei, which is largely unknown at this time. The location 
of the neutron dripline has implications for the nucle- 
osynthesis of heavy elements in the universe as well as 
for the physics of neutron stars. Radioactive beam facil- 
ities of the next generation are designed to establish the 
limits of neutron binding as one of their central goals. 
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We therefore want to discuss some of the straightfor- 
ward implications of our Woods-Saxon parameterization 
for the limits of nuclear binding. It is clear that a model 
like the Woods Saxon-potential does not include residual 
interactions and therefore a lot of interesting physics is 
out of reach for this approach. Nevertheless it can serve 
as the first order approximation in the determination of 
the drip lines and the basis for more detailed theoretical 
investigations. 

The calculated proton and neutron drip-lines are dis- 
played in Fig [11] and clearly show the effect of shell and 
sub-shell closures. Interestingly, the neutron drip-line 
shows a number of instances, e.g. the Ca isotopes, where 
the presence of magic numbers leads to regions, where 
an element has bound nuclides after unbound ones. This 
"meandering" of the neutron drip line does not occur for 
the proton drip line. 

Oxygen is the heaviest element, for which the location 
of the neutron drip line is known at the present time. 
Experiments with high sensitivity have established, that 
26 O is unbound [H, 63]. Our calculation, represented in 
Fig. QT]also shows 26 O to be unbound, but brings back 
28 as a nuclide bound by 0.21 MeV due to the N=20 
shell closure. 

The Tin isotopes are calculated to be neutron-bound 
up to 176 Sn. We display the calculated and experimen- 
tal neutron separation energies for the odd-neutron Sn 
isotopes in FigJTS] A surprising property of the calcu- 
lated energies is the almost constant, very low neutron 
separation energy of the Fermi orbital from N=93 up to 
the N=126 shell closure. This property is based on a 
large number of orbitals being situated very close to the 
Fermi energy- 176 Sn binds 34 neutrons within an energy 
interval of 0.84 MeV. For comparison, the same number 
of neutrons in the same orbitals is spread out over an 
energy interval of 3 MeV in the spectrum of 208 Pb. The 
large number of orbitals in a small energy interval will 
make it very difficult to predict the neutron drip-line in 
the Sn-isotopes by theoretical means and shows the need 
for experiments with very exotic nuclides. The situation 
of orbitals is analogous to the one in neutron rich Zir- 
conium isotopes with A>122, which was shown to lead 
to giant halos in relativistic Hartree-Boguliubov calcula- 
tions (r34| . This property should make the neutron-rich 
Sn nuclides a very interesting system to study new modes 
of collectivity based on the interaction with the contin- 
uum. 



VII. CONCLUSION 

We have established a new parameterization for the 
Woods-Saxon potential. Its six parameters are fitted to 



single-particle spectra around doubly magic nuclides and 
experimental charge radii. We achieved a high quality of 
description of the single-particle spectra observed around 
the doubly-magic nuclei between 16 and 208 Pb. The 
nuclear properties such as radii, shell evolutions, magic 
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Figure 12: Neutron separation energies for odd-N Tin iso- 
topes, calculated with the "Seminole" parameterization of the 
Woods-Saxon potential and experimental values from [HoT ] . 



numbers and driplines are all shown to be well described 
considering limitations of the Woods-Saxon approach. 
The isospin conservation by the nuclear part of interac- 
tion and proper treatment of the two-body particle-core 
kinematics are important for the success of the model. 
The parameter set is intended to enable the use of the 
Woods-Saxon potential as a basis for shell model calcu- 
lations and a pathway to connect the physics of bound 
and unbound nuclear states. 
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